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ABSTRACT 

We have investigated the short term linear relation between the rms variability and the flux 
in 1,961 observations of 9 black hole X-ray binaries. The rms-flux relation for the 1-10 Hz 
range is ubiquitously observed in any observation with good variability signal to noise (>3% 
1-10 Hz fractional rm s). This concurs with results from a previous study of Cygnus X-1 
jGleissner et al. 1 12004 . and extends detection of the rms-flux relation to a wider range of 
states. We find a strong dependence of the flux intercept of the rms-flux relation on source 
state; as the source transitions from the hard state into the hard intermediate state the intercept 
becomes strongly positive. We find little evidence for flux dependence of the broad-band 
noise within the PSD shape, excepting a small subset of observations from one object in an 
anomalous soft-state. We speculate that the ubiquitous linear rms-flux relation in the broad 
band noise of this sample, representing a range of different states and objects, indicates that 
its formation mechanism is an essential property of the luminous accretion flow around black 
holes. 

Key words: X-rays:general - X-rays:binaries 



1 INTRODUCTION 

Accreting compact objects such as black hole X-ray binary systems 
(BHBs) are known to display strong short-term variability, which is 
dependent on the accretion state. Systematic changes in the power 
spectra are observed as a source movesthrough different states over 
the course of an outburst (Belloni 20 IQ). The similarities between 
these changes for different sources suggest a common production 
mechanism, so investigation and comparison of the short-term tim- 
ing properties of different BHB systems is an important tool in un- 
derstanding accretion in these sources. 

One property which appears to be common to many accreting 
compact objects is a linear relationship between the rms variabil- 
ity of a source and the count rate. Initially discovered in observa- 
tions of one black hole and one neut ron star X-ray binary (Cy g X- 1 
and SAX 1808.4-3658 respectivelv: IUttlev & McHardvl200lh . this 
rms-flux relation has since been studied in a few other observations 
of BHBs, AGN, a neutron star X-ray binary and a ULX (see e.g. 
Uttley & McHardv"200l'; 'Uttlev' 2004: 'Gaskelll l2004l : lUttlev et al.l 
2005; Heil & Vaughan 2010; Heil et al. 201 1) and studied in detail 
for all Rossi X-ray Ti ming Explorer (RXTE ) observations of Cyg 
X-1 up to early 2003 jGleissner et ai]|2004l) . Optical observations 
of three BHBs in the l ow-hard state have also shown the rms-flux 
relation l lGandhill200t ). The range of compact objects from which 
this relationship has now been observed is striking and suggests 
that it could be a ubiquitous property of luminous accretion flows. 
Before this can be established it needs to be observed in a much 
wider range of sources and observations. 

The rms-flux relation in Cygnus X-I has been detected over 



timescales from seconds to years jGleissner et ai]|2004l) . This con- 
sistency and range of timescales is difficult to explain using models 
employing additions of independent shots; a connection is required 
between variability occurring over timescales of years and that hap- 
pening in sec onds. A propagati ng fluctuations model, such as that 
suggested bv iLvubarskiil ( Il997l) . where long-term variations from 
the outer reaches of the disc propagate inwards, coupling to short- 
term variations at smaller annuli, provides a natural explanation for 
the observed correlation dUttlev & McHardvll200lh . This modula- 
tion implies that the flux dist ribution from the source should be 
log-normal jUttlev et al.ll2005h . Log-normal flux distri butions have 
been observed from various BHBs and AGN (s ee e.g. lUttlev et all 
l2005l : lGaskellll2004l : ICiebels & Degrangell2009h . 

Although the short-term relation has been studied in many 
observations of Cyg X-1 and a selected few other BHBs, it has 
not been consistently observed in the full range of states exhib- 
ited by transient sources. The rms-flux relation may be ubiquitously 
present in all states, or alternatively the relation may not be seen in 
a particular state suggesting that there is something essentially dif- 
ferent in the nature of the accretion flow in between the states. In 
order to test the ubiquity of the rms-flux relation in BHBs we have 
utilised the full range of data from 9 BHBs in the RXTE archive. 



2 DATA ANALYSIS 

2.1 Lightcurve extraction 

We have utilised data from the RXTE archive for the 9 BHBs given 
in Table [T] taken by the Proportional Counter Array (PCA) until 
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Target 


Total obs. 


Obs. >3% CTfi-ac 


Good obs. 


Obs. IDs. 


(1) 


(2) 


(3) 


(4) 


(5) 


GX 339-4 


487 


376 


235 


P20056, P20181, P30165, P40031, P40105, P40108, P50117, P70109, P90704, P92035, 










P92052, P92085, P92428, P92704, P93409, P93702, P94331, P95409 (until 27th May 2010) 


V 1 '1 T 1 1 1 O 1 A Of\ 

Alh J 1 1 15+480 


48 


30 


26 


P50133, P50407 


GS 1354-64 


8 


7 


7 


P20431,P30401 


4U 1543-475 


97 


25 


12 


P70124. P70128. P70133 


XTEJ1550-564 


153 


106 


81 


P40401. P50134. P50135. P50137. P60428, P80135 


XTE J 1650-500 


177 


69 


39 


P60113, P60118, P70124 


GRO J 1655-40 


462 


140 


93 


P10261, P20402, P90058, P90428, P91404, P91702, P91704 


HI 743-322 


406 


297 


129 


P80135, P80137, P80144, P80146, P90058, P901 15, P90421, P91050, P91428, P92047, 










P93427, P94431,P95405 


XTEJ1859+226 


123 


55 


12 


P40122, P40124, P40440, P50401 


Total 


1961 


1105 


634 





Table 1. Note: The 1998 outburst of XTE J1550-564 is excluded from this data set due to the strong QPO; The 2004 outburst of GX 339-4 and 2005 outburst 
of XTE J 1 1 1 8+480 are mostly excluded as there are very few observations in the modes processed for this work. 



May 2010. The analysis was limited to those observations taken 
in either Binned, Single Bit or Event modes; observations with no 
data in any of these modes have been excluded from the data set. 
If observations were made with some or all of these modes in the 
energy range required, then Binned or Event mode data was used 
in preference to Single Bit mode to allow finer control over the en- 
ergy range selection. We have also excluded observations without 
a file in Standard 2 mode, as these data were used for the back- 
ground subtraction (this mainly applies to a few observations taken 
when there were problems with the hardware on RXTE). Proposal 
Identifiers with fewer than six observations are also excluded from 
the data set. These objects have been selected as they are all well 
known black hole X-ray binaries, which are bright in outburst and 
have a number of observations within the RXTE archive. This is 
not a complete sample of low mass BHBs in the RXTE archive, 
but is chosen to be representative of these sources in general. All 
lightcurves have been extracted and then binned to around 3 ms 
(depending slightly on the time resolution of the data mode). In 
order to adjust for the differing combinations of detectors used in 
observations taken with the PCA all results presented here are nor- 
malised to one PCU (proportional counter unit). All of these obser- 
vations are taken when the object is in outburst; in quiescence the 
count rates are generally too low for detailed timing analysis. Con- 
sequently the observations are not evenly sampled in time but do 
represent a range of states for each source. Due to problems caused 
by a strong quasi-periodic oscillation (QPO - discussed in Section 
[33] l most observations from t he 1998 oiitburst of XTE J1550-564 
are excluded from the sample jHeil et al]|201lh . 



In general the analysis was performed on data extracted from 
the energy range closest to 2-13 keV, but this differs slightly be- 
tween observations due to the different data modes used and vari- 
ations in detector operation. As Standard 2 channel is thought 
to introduce artificia l variability into the lightcurves (see e.g. 
iGleissner et"aLll2004l) this channel was excluded. The lightcurves 
were then background corrected using simulated background 
lightcurves produced by the ftool PCABACKEST. We applied the 
bright source model to all observations, despite some having count 
rates below the 40 ct/s/PCU threshold where the "faint source" 
model is usually preferred. Our choice ensures some consistency 
in the analysis of observations on either side of the threshold, but 
experiments using both background models confirm that the choice 



made little difference to the overall results Observations with a 
count rate < 10 ct/s/PCU are excluded from the sample. 

2.2 Calculation of the rms-flux relation 

The resultant background corrected lightcurves were then split 
into contiguous segments of 3 s in length and the mean count 
rate ((F)) was found for each segment. Periodograms were then 
measured for the segmen ts using absolute normalisation (see e.g. 
lUttlev & McHardvll200l[) . in this normalisation the integrated area 
underneath the power spectrum is equal to the variance. The 
stochastic nature of the lightcurve results in scatter in the peri- 
odogram, the measured power values are scattered around the true 
value acc ording to a dist ribution with two degrees of freedom 
(see e.g. Ivan der Kli j[T989h . In order to reduce th s effect and 
estimate the real value of the power in each frequency bin, the 
periodograms were averaged into flux bins. The number of peri- 
odograms per flux bin was dependent on the number of segments 
in the individual observations but ranged between 25 - 100, ensur- 
ing each observation contained a minimum of 10 flux bins. In order 
to find the rms, we first found the total variance in the required fre- 
quency range, given by summing the power over vi - V2, and then 
multiplied by the frequency resolution (A/) i.e. 



(1) 



Where 5 is variance, (yj) is the power in a frequency bin with 
the bar indicating that this is averaged, and W is the number of peri- 
odogram points in the frequency range. In this paper the frequency 
range was chosen to focus on the region of the power spectra which 
shows power consistently in every state, whilst also displaying 
some of the most interesting features of the power spectra (spec- 
tral breaks and QP Os) - 1-10 Hz. It is the convention with rms-flux 
relations (see e.g. [Utflev & McHardvll200ll ; IGleissner et~al]|2004l ; 



lUttlev et alj|2005l : 



Gandhill2009^ to sample relatively restricted (~ 



1 decade) frequency ranges for a number of reasons: Firstly, a lower 
limit of 1 Hz allows us to use shorter segments of lightcurve (in this 
case 3s), allowing us to measure the rms-flux relation on shorter 
timescales and to retrieve more flux bins from shorter observations. 
This means that a larger flux range is measured and the shape of 
the rms-flux relation is better defined. There is also a consistently 
high level of power at 1 Hz, particularly in the harder states (see 
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Figure 1. Examples of noise-subtracted power spectra for the upper and lower flux quaililes and their corresponding rms-flux relations. These are normalised 
to the rms normalisation used to calculate the rms-flux relation, for this reason the flux dependence is still clearly visible in the power spectra. Low flux quartile 
power spectra are green, high flux are purple. Overleaf: A: XTE J1550-564 (obs.id. 80135-01-04-00; Hard State) B: GX339-4 (obs. id. 95409-01-12-01; Hard 
State), CXTEJ1550-564 (obs. id. 30188-06-01-02; HIMS)Z): GX 339-4 (obs.id. 70109-04-01-01; mUS) Above: E: GRO J1655-40 (obs. id. 91702-01-63-00; 
Anomalous Soft State), F GRO J1655-40 (obs. id. 91702-01-62-00; Anomalous Soft State) 



e.g. iGierlinski et alj|2008h. which again means we sample a wider 
ranee of fluxes. IGierlinski et all ( |2008|) show that the broad band 
noise displays little evolution within this frequency band, meaning 
that the stationarity in PSD shape required to measure the short- 
term rms-flux relation is more likely to be observed - this excludes 
certain observations with QPOs within this frequency range, which 
may display different behaviour to the unde rlying noise (see e.g. 
iHeil et aflllOllI : [Rodriguez & VarnierelllOllI) . Figured shows that 
although there is often some power above ~10 Hz, in general this 
is insignificant when compared to that observed within the chosen 
1-10 Hz range. 

As the power in each frequency bin contains a contribution 
from the noise this needs to be estimated and removed so that the 
rms can be found & = - Pf^W^f, where is the power den- 
sity for the Poisson noise component of the variance and is given 
by Pn = 2(x -I- B), X is the source count rate and B is the back- 
ground. The flux binning allowed for an accurate calculation of the 
errors on the binned & va lues, these were found in a more general 
manner to that applied in lOleissner et al See appendix A 

for further discussion of this error. 



2.3 Fitting the linear relation 

Following lUttlev & McHardvl ( 1200 ih we have assumed that the 
rms-flux relation is linear and described by 



& = k(F- C,-) 



(2) 



Where the gradient of the relation is given by k, and C, is the inter- 
cept on the X-axis. Non-zero inte rcepts have been widely observed 
previously. iGleissner et al. I( l2004l) found both positive and negative 
X-axis intercepts in the rms-flux observations from Cygnus X-1, 
possibly indicating the presence of two variable components within 
the light-curve. 



2.4 Testing for "good" rms-flux relations 

In order to evaluate the rms-flux relation in each i ndividual obser- 
vation we employed a method similar to that used in lOleissner et al.l 
( |2004t) . Test ing for correlation using Kendall's rank correlation co- 
efficient, T jPress et al]|l992l) . in addition to the test to eval- 
uate the fit. Kendall's tau tests for correlation which the bulk of 
points follow without being heavily aff^ected by outliers. Due to the 
lower levels of variability and count rates in data from the BHBs 
compared to observations of Cygnus X-1, less stringent limits 
were u sed to define "good" data than those used in lOleissner et al.l 
| |2004|) . Whilst their minimum threshold was set at t„,„ = 0.9 we 
have used t,„,„ = 0.5. Observations with this value of t still display 
a positive linear correlation but greater scatter is allowed around the 
trend. The assumption was made that if r was significant at a level 
of IcTj the data were sufficiently strongly correlated. In addition to 
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this we require the fit to the linear model to be accepted using a < 
0.003 (i.e. not rejected at Str^z). 

In general the data were found to agree with the assumed lin- 
ear model, with the majority of discrepant observations due to data 
at low count rates and/or with low fractional rms. Assuming a back- 
ground rate of approximately 7 ct/s/PCU and a count rate of 100 
ct/s/PCU the fractional rms of the Poisson noise level is 2.14%. Al- 
though this level drops rapidly as the count rate increases, a signif- 
icant amount of power above the noise is required to satisfactorily 
measure the rms. Some observations with count rates below ~ 500 
ct/s/PCU and a fractional rms of ~3% (where the noise level drops 
to 0.4%) may only just have enough power for the source vari- 
ance to be measured accurately in the observed frequency band. For 
this reason a limit of 3% fractional rms was imposed, observations 
where the fractional rms is less than this within the frequency band 
were discarded. Most observations with count rates lower than 100 
ct/s/PCU are typically observed in a hard state where the fractiona l 
rms is higher than 10%c (see e.g. iRemillard & McClintoc3l2006l) 
and so are less strongly affected by this limit. However, in a stan- 
dard soft state the variability often drops drop below the required 
3% fractional rms level, so even though the count rate is high few 
observations of sources in this state are included in the following 
analysis. 

Many of the other rejected observations had low count rates, 
resulting in flux binned PSDs with low signal to noise and negative 
Poisson noise-subtracted variances. Observations containing nega- 
tive variance bins were therefore excluded from the sample. Short 
observations also limit the signal to noise for the flux binned PSDs 
as the number of rms points per flux bin are low. The length of 
good time required for a detected rms-flux relation varied accord- 
ing to the count rate and the level of intrinsic variability within the 
observation, but more than 1 ks was generally required. 



3 RESULTS 

3.1 Ubiquity of the rms-flux relation 

1,961 observations of the sources listed in table [T]have been tested 
for a simple linear rms-flux relation. Of these observations 1105 
displayed the 3% fractional rms required for the rms-flux relation to 
be measured. A "good" rms-flux relation as defined by the criteria 
set out in section l2!4l was observed in 634 of these observations. 

Observations where a positive rms-flux correlation was not 
found were inspected individually. Non-stationarity of the power 
spectrum within an observation can cause a 'wavy' rather than a 
linear r ms-flux relation to be o bserved, such as those shown in fig- 
ure 8 of iGleissner et al.l ( |2004|) . To some extent a further example 
is shown in Figure[T]panel E, there is a lower flux segment of time 
series within this observation where the rms-flux relation is slightly 
flatter, for this reason the lowest flux bin has a slightly higher rms 
when compared to the other bins. The flux-binned PSDs of long, 
bright observations which deviated significantly from a positive lin- 
ear rms-flux correlation were inspected for any shape change over 
the course of the observation. In some observations plotting the in- 
dividual flux binned power spectra for each point in the rms-flux 
relation, revealed QPOs shifting with flux in either frequency or 
amplitude over the course o f the observation, similar to those dis- 
cussed in iHeiletall ( 1201 ih and iRodriguez & Vamiir3 ( 1201 ll) . In 
other cases it was apparent that the shape of the whole power spec- 
trum was evolving with flux. No examples of non-linear or neg- 
atively correlated rms-flux relations were found within the sample 
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Figure 2. Results for all good mis-flux relations. The colours correspond 
to the state of the source and states are detined as follows: LHS - Low 
hard state; HIMS - Hard Intermediate State; SIMS - Soft InteiTnediate State; 
Anom SS - Anomalous Soft State; HIMS/Anom SS - observations with char- 
acteristics of both the HIMS and the Anomalous Soft State which are 
uniquely observed in H 1743-322. The contours shown are lines of frac- 
tional i-ms at 30%, 10% and 3%.. 

that could not be explained by one of the above reasons. Where suf- 
ficient data are available and the variability is stationary, the rms- 
flux relation appears to be ubiquitous in the broad-band noise from 
BHBs. 

Figure[T]illustrates observations of the rms-flux relation across 
a range of different states and sources, from the low-hard to an 
anomalous soft state. 

Figure [2] shows the gradients against the intercepts for all 
sources within the sample. As the flux-intercept (Cj,) is dependent 
upon source count rate the values plotted are shown as a fraction of 
the mean count rate over the observation (the fractional intercept: 
C.J{F) = Cfrac), this allows the parameters from all sources in all 
states to be compared directly. The "banana" shaped paths within 
this plot are explained as lines of fractional rms if we re-write the 
gradient as 



In this case the value of the gradient may be interpreted as the frac- 
tional rms of the observation after removal of the intercept on the 
X-axis, if C, = then k - {a-fmc) (where {&fr„c) = (o">/(F) and 
is fractional rms). This illustrates the importance of the intercept in 
interpreting these results. Observations with the same overall frac- 
tional rms and count rate but different intercepts (C^) will show 
rms-flux relations that pivot around a common point, meaning ei- 
ther steep or shallow gradients may arise from observations with 
the same mean rms and flux. When k is plotted against Cfmc then 
the changes in C,- and (F) work to produce the observed curves: 
— » oo as Cfrac — * 1 and A; — » as Cfrac The lower and 

further to the right the curve on Figure (2] the lower the fractional 
rms of the observations. As f ractional rms is well known to b e re- 
lated to source state (see e.g. iMcClintock & Remillardl2006h then 
the position of points on the k vs. Cfrac diagram provides a link to 
the energy spectral properties of the observation. 

Whilst the paths on this plot are therefore just contours of dif- 
ferent fractional rms values, it clearly displays the full range of 
intercepts and gradients seen within the sample. It is also appar- 
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longer dmescales there is not the equivalent stationarity in power 
spectral shapes observed within a single observation, therefore the 
long term relationship in these states is not representative of that 
observed on short timescales. In contrast evolution of the power 
spectra within the hard state occurs mainly below the 1- 10 Hz range 
and although the shape does not remain completely stationary it is 
more suitable for comparison with the short term rms-flux relation. 

In Cygnus X-1 separate linear rms-flux relations were ob- 
served in the long-term relations in both the hard and softer states. 
Figure [3] shows the mean rms and flux values for all good rms- 
flux observations within the sample. There are two clear paths, the 
upper one is that followed by sources in the hard state, the lower 
shows those in anomalous soft states. Intermediate states fall be- 
tween these two lines. The lower edge of this plot is artificially 
introduced by the requirement that there be at least 3% fractional 
rms within the 1-10 Hz frequency band (marked by the dotted line), 
so there is a wider range in the lower fractional rms soft state ob- 
servations than appears here. However, the relations for individual 
outbursts appear to be remarkably linear over a wide range of count 
rates, and there is distinct similarity in those observed between 
different sources. The agreement in these pat hs over repeated out- 
bursts of GX 339-4 has been remarked upon in lMunoz-Darias et al.l 

If) and we return to this point in section [33] This plot also 
clearly shows the sharp drop in fractional rms over a few obser- 
vations within the transition from the HIMS through the SIMS to 
soft anomalous states, seen in the points which lie in between the 
hard (upper) and soft (lower) state rms-flux relations. The lack of 
intermediate points at high count rates occurs because these obser- 
vations are taken during bright flares where no state changes were 
observed. Further differences between the relations over the course 
of the outbursts are discussed in the following section. 



Figure 3. Top. Long term rms-flux results for all good rms-flux relations in 
the 1-10 Hz range. Bottom As Top but concentrating on lower values. The 
dotted line plots the lower limit on detection for the rms-flux relation - 3% 
fractional rms. The colours correspond to those used for Figure|2] 



ent that is less well constrained at negative values than posi- 
tive ones, this is partly an artefact of the fitting, by nature nega- 
tive intercepts are further from the actual observed rms-flux values 
than positive ones and therefore the point where the fitted line in- 
tercepts the X-axis is less certain. The range of intercepts and gradi- 
ents observed within the sa mple is striking, partic ularly when com- 
pared to those observed bv lGleissner et al.l ( |2004|) from Cygnus X- 
I. Whereas their figure 5 indicated only a few observations with 
negative flux-intercepts in the softer states and none in hard states, 
we observe a wide distribution of positive and negative values 
over a range of fractional rms values. Section [331 explores the po- 
sitions of points in Figure|2]in relation to source state. 



3.2 Long term rms-flux relation 

In a similar manner to iOleissner et alj j2004h the long-term rms- 
flux relation is plotted for all of the sources in the sample by tak- 
ing the mean rms and flux values for each observation in the 1-10 
Hz band. The rms-flux relation measured within each observation 
relies on the PSD remaining the same shape in the 1-10 Hz fre- 
quency band. The power spectra in intermediate states can show 
dramatic changes in between observations, with QPOs and broad 
band components changing in both frequency and amplitude. On 



3.3 Evolution of rms-flux relationship with BHB state 

lOleissner et"ai] ( |2004|) noted a distinct difference between the be- 
haviour of k and C, in different states of Cygnus X- 1 . In the softer 
states k was generally found to be lower and Cx had a much wider 
range of values than those observed in the hard state. 

Figure [4] shows the evolution of the hardness ratio, fractional 
rms, gradient {k) and intercept iCfrac) for sources with the largest 
number of good observations. The flux-intercept is again shown as 
a fraction of the mean count rate, allowing for comparison between 
different sources and different states of the same source. All four 
parameters plotted are dimensionless. 

The state of a source is usually identified based on a com- 
bination of its luminosity, timing properties (the power spectrum) 
and energy spectrum. The most important single diagram cur- 
rently in use for diagnosing and tracking the evolution of states 
is the hardness intensity diagram (HID). The outbursts of tran- 
sient black holes, and some neutron star and cataclysmic vari- 
able syst ems tend to follow a similar "q" shape path around this 
diagram (iHoman et al.l 1200 ll : iBelloni et al] 120051 : iBellonil |20IO| : 
iKording et al.ll2007l) beginning in a low/hard state (HS), evolving 
through hard then soft intermediate states (HIMS/SIMS) to the soft 
(thermal dominated) state, before changing back to the hard state 
via lower luminosity intermediate states. The HID itself is a pow- 
erful tool based on only the source intensity (in a particular X- 
ray band) and a very crude indicator of the energy spectral shape, 
w ith no explicit d ependence on tim ing p roperties. But a s discussed 
in lMcClintock & Remillarj | |2006|) and iBellonil ( |2010|) the timing 
properties of these objects evolve as clearly, if not more clearly 
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Figure 4. Hardness ratio, fractional rms, Gradient and Intercept vs. Count rate for four sources. Points below 3% fractional rms are excluded from the analysis. 
Colours correspond to those used in Figure[2] 
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than the energy spectral properties do. For more details see lBellonil 
.Dunn et al.. (.20m) and references therein. 
It is already well-known that the fractional rms (in the usual 
2-20 keV X-ray band) i s lower in the soft than hard state. 
iMufioz-Darias et ^ ( 1201 ih presented a new plot which demon- 
strated how the pattern described by the HID is repeated by 
the rms. In the second panels of Figure |4] we show an altered 
version of this new diagram - the rms intensity diagram (RID) 
- the natural variability counterpart to the HID, which shows 
the intensity against fract ional rms amplitude as discussed in 
iMufioz-Darias et ^ | |201 ih . There is a remarkable similarity be- 
tween the HID and RID. The RID is of course very similar to the 
long-term rms-flux plot in Figure (3] indeed it is simply a geomet- 
ric tr ansformat ion of that p lot and contains the same information. 
iMufi oz-Darias et alj ( l201lh clearly demonstrate that the position on 
the RID is related directly to the source state in GX 339-4. This 
demonstrates that variability is every bit as strong an indicator (or 
diagnostic) of outburst evolution and states as is the energy spec- 
trum. 

For GX 339-4, XTE J1550-564 and XTE J1650-500 the clas- 
sic 'q' shape of the hysteresis is replicated by the fractional inter- 
cept over the course of the outburst. The intercept appears close to 
zero in the hard state, but becomes strongly positive as the source 
transitions into the HIMS before moving back again. This can be 
seen in Figure [2] the red and blue points (HS and HIMS respec- 
tively) occupy a similar area within the plot (which is expected 
as one evolves from the other) but the highest values of Cfi.„c oc- 
cur when the source is in the HIMS. This behaviour is seen in at 
least 3 repeated outbursts of GX 339-4 which all follow remarkably 
similar paths in both the hardness and rms - intensity diagrams. 
In contrast the high flux observations of GRO J 1655-40 show a 
wide range of intercepts including ones which are strongly nega- 
tive. In this case the power spectra indicate that the source is in an 
anomalous soft state teellonfcOld) , suggesting that the behaviour 
of both the intercept and gradient is strongly state related (these 
observations correspond to some of the green points in Figure O. 
The observations of GRO J 1655-40 in the hard state and HIMS still 
display similar strong positive intercepts observed in other sources. 
This can be observed in the lower flux (< 500 ct/s/PCU) observa- 
tions in Figure |4l The appearance of intercepts close to zero in the 
hard state rising to strong positive values in the HIMS and then dis- 
playing a wide range of values, including strong negatives, in the 
soft states is replicated for other sources in the sample (i.e. H1743- 
322). 

The state dependence of the intercepts, explains the wide 
range of values in Figure |2] at low fractional rms, and almost bi- 
modal spread between those with high and low k values. Observa- 
tions with Cfrac > 0.5 and k > 0.1 are from a hard-intermediate 
state (blue), whereas when Cf,-„c < 0.5 and k < 0.1 the source is in 
an anomalous soft state (green). Observations in the SIMS have a 
range of values but this is likely to be becau se they are strongly in- 
fluenc ed by the behaviour of the QPOs (see lRodriguez & Vamierd 
I2OI ih . Figure[2]demonstrates that both k and C/iac are not only de- 
pendent on the fractional rms within the 1-10 Hz range, but also 
upon the state of the source. 

3.4 Flux dependence of the PSD shape 

In order to investigate any possible flux dependence of the PSD, 
power spectra were measured for 100 s long segments of lightcurve 
and grouped into four flux bins. The averaged PSDs from the high- 
est flux bin were then compared with that from the lowest. Figure 



[T]shows the results for the six example observations. We tested the 
constancy of the PSD shape by fitting the ratio of high and low flux 
PSDs with a constant, using the statistic to evaluate the quality 
of the fit. The ratio was fitted up to 20.0 Hz, above this point some 
observations are dominated by Poisson noise (see e.g Figure[Tf). If 
the fit was judged to be poor (> 3(r^2 ) then the observations were 
subject to further investigation. 

Little evidence for flux dependence was found over the ob- 
served frequency range in most of the observations: in total only 
20 of the good observations were identified as having obvious 
frequency-dependent differences between the PSDs from the up- 
per and lower flux quartiles. Different binning was tested in order 
to eliminate possible spurious results caused by sharp features (e.g. 
QPOs) appearing shifted in frequency between the two power spec- 
tra. Fifteen of these observations were from the 2005 outburst of 
GRO J 1655-40 whilst it was in the anomalous soft state. The power 
spectra all show strong low frequency noise with a high frequency 
broad band component, an example of which is shown in Figure[l^. 
The two flux binned power spectra indicate that the peaked band- 
limited component shows a greater change in rms with flux than 
the power law. The observations with clear flux dependence are not 
consecutive and the appearance of the power spectrum changes be- 
tween this and the typical soft state shape. The power spectra in 
Figures [T^ and f were taken almost exactly a day apart and show 
distinctly different behaviour, although the fractional rms in the 1- 
10 Hz band remains similar (3.7 % vs. 3.5 %). Figure[Tf is typical 
of that observed in the soft state, only with a higher fractional rms, 
whereas Figure[T^ is more "anomalous". 

The upper and lower flux binned power spectra in Figure [1^ 
have been separately fitted using xspec to test if the change in 
power spectral shape can be mainly explained by the strength of the 
peaked noise components. The model used a bending power, plus 
two Lorentzians for the high frequency broad band component and 
a separate narrow Lorentzian to model the very low frequency com- 
ponent, which can be seen in the lowest frequency bin of Figure[T^. 
This model gave a fit to both the high and low flux power spectra 
with good X' statistics (x' = 35 and 47 respectively with 38 dof). 
The models used were broadly similar for both of the spectra. The 
main difference was in the first broad Lorentzian at high frequen- 
cies which was wider in the high flux spectrum, (2 = 1-9 compared 
to 1.0 in the low flux spectrum). Other than this the only difference 
between the two models was in the normalisation of the compo- 
nents. The drop in the rms for the power law component between 
the high and low flux power spectra was 21%, in contrast the rms 
of both the combined Lorentzians for the high frequency peaked 
component and the Lorentzian for the lowest frequency component 
dropped by 67%. This suggests the existence of two separate com- 
ponents combining to produce the observed light curve, one form- 
ing the peaked noise and another the underlying broad-band com- 
ponent. Both of these components show a positive linear rms-flux 
relation, but have different flux-axis intercepts. This might help ex- 
plain the well constrained negative intercepts observed in the rms- 
flux relations in some power spectra from this state (see Section 
[33}. 

The other 5 observations which did not meet the criteria for 
inclusion (i.e. did not display a consistent flux dependence over 
the whole frequency range) were generally found to contain com- 
ponents shifting in frequency over the course of the observation, 
sometimes outside the 1-10 Hz range. These shifts, which have in 
these cases still produced strong, positively correlated, linear rms- 
flux relations, demonstrate the importance of careful analysis of the 
power spectra when measuring the rms-flux relation. Small differ- 
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ences in the shape of the PSD, whilst not always drawing it away 
from linearity, may artificially alter the gradient and thus the inter- 
cept of the observed rms-flux relation. 



4 DISCUSSION 
4.1 Summary 

A large sample of observations from nine transient black hole X- 
ray binaries have been tested for a linear rms-flux relation. Analysis 
of the sample revealed the following results: 

• In every observation with sufficient flux and intrinsic variabil- 
ity to measure the rms, and a stationary PSD, a positive linear rms- 
flux correlation was observed in the 1-10 Hz range. 

• All sources and states show long-timescale rms-flux rela- 
tions (for observations with little alteration to power spectral shape 
within the 1-10 Hz band), similar to that already seen in Cygnus 
X-1. These bear clos e relation to those observed in Cygnus X-1 
dOleissner et al.ll2004l) and show remarkable similarities both in re- 
peating outbursts of the same source, and in the behaviour of dif- 
ferent sources (i.e. lMufioz-Darias et alj|2oTll) . 

• There appears to be a link between k, Cx and source state. A 
wide range is observed in the flux intercept (C^). Hard states have 
X-axis flux intercepts close to zero, strong positive intercepts are 
observed in the HIMS. Strong negative intercept seem primarily 
isolated to the soft states and it is for observations in these states 
that the greatest range in C j is observed. 

• There is little evidence for strong flux dependence in the PSD 
shapes apart from in some observations of an anomalous soft state 
of GRO J1655-40. 



4.2 The ubiquity of the rms-flux relation 

The apparent ubiquity of the rms-flux relation in the broad-band 
noise within observations of these black holes implies that the cou- 
pling of variability over all timescales is common to the accretion 

process in all states. 

Models similar to that explored in iLvubarskiil l ll997h where 
perturbations in the accretion rate propagate inwards through the 
accretion disc could provide an explanation for this behaviour. 
Variations on long timescales originating at larger radii in the 
accretion disc propagate inwards, coupling with faster perturba- 
tions at smaller radii, therefore creating a continuous relation be- 
tween the mean accretion rate and amplitude of variations on all 
time scales, variati ons r eaching inner regions mo dulate the X-rays. 
lUttlev et af] ( l201ll) and iRevnivtsev et alj j201lh have studied ob- 
servational evi dence for the propa gation of variability through the 
accretion disc. lUttlev et al.U201lh use data from GX 339-4 taken 
with XMM-Newton which reaches softer energies than observa- 
tions with RXTE, in order to test for time lags between different en- 
ergy bands. They find clear evidence that soft disc black body emis- 
sion leads hard emission over long timescales (> 1 s). The length 
of the observed lags suggest that they are generated by the inward 
prop agation of accretion fluct uations through the X-ray emitting re- 
gion. IRevnivtsev et alj ( 1201 ih observe some optical emission lead- 
ing the X-rays in simultaneous observations of three intermediate 
polars, the lead-time (~ 7 s) agrees closely with that expected for 
the travel time of captured material onto the white dwarf. These 
results are part of an increasing array of evidence for the propagat- 
ing fluctuations model, and suggest that a linear rms-flux relation 
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Figure 6. HID (left) and RID (right) for GX 339-4. Observations shown are 
those from the 2007 outburst and beginning of that from 2010 which have 
count rates less that 100 ct s"'. Points are colour coded according to MJD 
from green to red. 

should be expected in the broad band noise within observations of 
a whole range of accreting compact objects. 

4.3 Interpreting state dependence of the intercepts 

Plotting k and Cf,c,c against count rate reveals clear state de- 
pendence, when compared to both the HID and RID (see Fig- 
urelUl. The RID shown he re is similar to th o se firs t presented in 
iMufioz-Darias et alj bOUl) and iMotta et all l |201l|) , in this case 
however the x-axis is (&frac), not simply (&). Figure|5]shows all ob- 
servations of GX 339-4 processed for this work with those used in 
the analysis marked in red. Lines o f constant fractional intercep t are 
diagonal in the plots presented bv lMuiloz-Darias et al.l ilOl ih and 
vertical here. The RID and the HID do not have identical shapes (al- 
though the mapping between the two is very close). Some features 
are not present in both plots - for example, the RIDs for both GX 
339-4 and XTE 11550-564 show separate paths in the RID which 
ar e not apparent in the HID. I n the case of GX 339-4, as described 
in |m ufioz-Darias et 311(1201 l|), the start of the outburst and return 
to quiescence are along separate paths in the RID (see Figure [6l( 
although they appear to merge as the source fades. This behaviour 
is observed over repeated outbursts but occurs at very similar val- 
ues of the hardness ratio and is therefore hidden in the HID. There 
is also a third path visible in Figure |6l followed when the source 
flux rises from low flux levels at the end of the 2007 outburst but 
then returns to quiescent levels. Similarly in XTE J1550-564 the 
return to the hard state following an outburst where state changes 
have occurred, is along a different route in the RID to that observed 
for outbursts when the source remains in the hard state. In the three 
hard outbursts XTE J 1550-564 rises and decays along very simi- 
lar paths of fractional rms. Again these tracks are much harder to 
discern on the HID. 

Conversely, there are features visible in the HID which are not 
as clear in the RID. For instance there are many good observations 
of GRO 11655-500 taken when the source was in an anomalous 
soft state; two separate types of power spectra are visible whilst the 
source is in this state (see e.g. [T^ and f). Those which display a 
typical soft state shape (i.e.[TJ) have a softer energy spectrum than 



© 2002 RAS, MNRAS OOO.fTlflll 



10 



n — I — I I I 



"T 1 1 1 — I — I I I 



n 1 1 — r 



1000 



■-H-::: 



*'f ' ■■ " 



3 

o 

a. 



A 
b. 
V 



100 - 



10 



^± .■ 

+ 



-+ + 



+ 



.•■.'+••■ 



_l I I I I I I I 



^1 



r 



_l I I L 



1 1 1 1 I I 



"T 1 1 1 1 — I I I 



I — ^ 



HH 



PHI-+fc|l..r^......-.+. 



.1 1 1 



I — I 1^ ■ ■* 



I'llll 



..-■■I - . -.. -•'< .. i:7T fp:. 



1 

I— — I 



I — — I 
I 




J I I I I I 



_] I I I I I I I 



0.01 0.10 

Hardness ((9-20keV)/(4-6keV)) 



0.01 



0.10 



Figure 5. HID (left) and RID (right) for all observations of GX 339-4. The "good" observations are shown with red symbols, and black symbols indicate the 
data (mostly soft states) with rms below our threshold for rms-flux analysis. 



the observations which are more intermediate (i.e.[T^), however the 
fractional rms in the 1-10 Hz range is very similar for both and 
this difference can only be observed in the HID. The use of both 
of these diagrams in conjunction allows for a strong identification 
of source state without the need to study details of the energy or 
power spectra individually. 

The well constrained positive and negative flux inter- 
cepts found for the rms-flux relation have been explained 
in prev ious work as secondary components wi t hin th e light 
curve. Juttiev & McHardyl 1200 ll: iGleissner et af] l2004h . Both 
iGleissner et all 1 2004) and Vaughan et all 1 201 ih find e vidence for 
energy dependence in the offset. Ivaughan et show that 

in observations of the AGN NGC 4051 the spectrum of the flux 
offset closely follows the soft component in the energy spectra, 
disappearing entirely at hard energies (> 2keV). This could sug- 
gest that the strong positive intercepts observed in HIMS within 
the sample are due to a soft non-varying component contributing 
with increasing strength as the source transitions into a softer state. 
However this does not explain the wide range of behaviour visi- 
ble in the anomalous soft states. Figure |4ji shows values between 



-0.98 < Cfrac < 0.5, although the gradients remain low (k < 0.1), 
any additional components would have to be varying in strength 
and power quite considerably to create this range. The lack of flux 
dependence within the PSD shapes for most observations within 
the sample is also important when considering the origin of any 
additional components. It suggests that the power spectra of the 
components must have similar shapes. The outlying cases (i.e Fig- 
ure [5) may be caused by either the presence of two components 
contributing on different timescales, or one additional component 
contributing to a small frequency range and reacting differently to a 
change in the flux. Loosely classifying the component's frequency 
ranges into 0.1-1.0 Hz and 1.0-20.0 Hz for the case in Figure[T^ the 
observed intercepts are very different. In the former range the frac- 
tional intercept is C^,„j.„j = -6.61 + 1.4, in the latter Cfrar|„_,„Q 
= 0.1+0.01 in both cases the likelihood of the intercept actually 
being zero is excluded at the 99% level. The different frequency 
ranges observed may also be indicative of the components origi- 
nating from different locations, as it suggests that they dominate 
over different timescales. This wildly differing behaviour demon- 
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strates that additional components within the lightcurve are likely 
to be present in at least some cases. 

This work demonstrates that the rms-flux relation is a ubiqui- 
tous property of the broad band noise in light curves from accreting 
black holes and that the behaviour of this relation is clearly de- 
pendent on the state of the source. Further work on other compact 
objects, i.e. neutron stars and cataclysmic variable stars, AGN and 
ULXs, could reveal that this is in fact a ubiquitous property of the 
luminous accretion flow in all compact objects. 
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APPENDIX A: ERROR CALCULATION FOR RMS-FLUX 
POINTS 

Here we derive the formula used to estimate the standard deviation 
of the estimates of the rms used in the rms-flux an alysis. We have 
calcu lated them in a more general manner to that of jCleissner et al.1 
l2004h . Each time series x,, sampled at a rate Af, is broken into seg- 
ments of length Wseg and for each segment a periodogram is com- 
puted, and the mean count rate estimated, (x'>. The periodogram 
values are calculated at Fourier frequencies fj = j/N^^g^t (for 
j = 1, 2, . . . , A'seg/2), and the periodogram value at the 7th fre- 
quency is y'j for segment /. The frequency resolution of such a peri- 
odogram is defined by the length of the segments: A/ = l/Nxu^t- 
The power spectrum at a given flux level is estimated by ar- 
ranging the periodograms in order of the corresponding average 
count rate (x'>, and averaging them in groups of M periodograms. 
For each count rate bin we calculate the average periodogram 

1 " 

It is well known that, in the absence of sampling distortions, 
the periodogram ordinates are distributed about the true spectrum 
following a chi square distribution (see Ijenkins & Watra | l96S ; 
lGroth"l975'; 'Priestley" 198 1'; 'Leahy et al." 1 9831: Ivan de r Klis"l989 ; 
(Press et al. 1992; Percival & Walden 1993: Timmer & K onig 1995; 
lBloomfieldl2000l ; IChatfielj|2003l) . Specifically, y. ~ Pjxl/2 where 
Pj = P(fj) is the true power density at frequency fj and xl is a 
random variable with a chi square distribution (with two degrees of 
freedom). 

The expectation and variance of the chi square random vari- 
able are Elxl] = 2 and yixl^ - 4- Using the rules for manipulating 
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expectations and variances (e.g. lEadie et al.ll97ll p23-24) we have 
E[y'j] = Pj and V[yp = P] and so 



1 



1=1 



(A2) 



(A3) 



where we have assumed the process is stationary, such that E[y'^] = 
Pji?, the same for all segments ; (within the given count rate bin), 
and the individual y'j values are distributed independently for all i. 
The above results show that the averaged periodogram is an unbi- 
ased estimator of the true spectrum P{fj) and its fractional "stan- 
dard error" is M 

An estimate of the variability power, 5, contained within a 
specific frequency range, say fj to // + {W - 1)A/, is given by 
summing the powers in the averaged periodogram {yj) over the ap- 
propriate frequencies: 



(A4) 



Given the above information it is straightforward to calculate the 
expectation and variance of the power estimate S : 



E[S] = A/2e[(>v)] = A/J]/'„ 



(A5) 



in all frequency bins, i.e. locally white noise. When this is not the 
case, e.g. when most of the power comes from a few dominant fre- 
quencies (strong power laws or peaked QPOs), both the power esti- 
mate and the variance on the power estimate are dominated by the 
frequencies that contribute the most power. One way to see this is to 
consider adding mor e frequencies th at contribute negligible power 
to the sum. Using the lGleissner et aL formula the relative er- 

ror on the total power is (MW)"''^ (from their equation 2), whereas 
the above equations predict this as (2 Pj/Af)''^/ 2 Pj (combining 
equations IA5l and lA6t . The former decreases as W increases even 
if the additional frequencies contribute no power, whereas the latter 
formulation leaves the relative error unchanged, as expected if the 
additio nal data contribute nothing. In general the lOleissner et all 
formula will underestimate the errors, i.e. produce confi- 
dence intervals that are smaller than the expected scatter in the rms 
data. 



V[5] = A/;^V[(v,)] = ^2/'J 



(A6) 



where we have assumed that the periodogram estimates, and their 
averages yj, are independently distributed for all j, which is asymp- 
totically true (as N^^s. ^ °°)- These results show S is an unbiased 
and consistent estimator of the integrated power spectrum over the 
required frequency range. 

The estimated power 5 includes a contribution from the Pois- 
son noise which is expected to contributed a flat spectrum with den- 
sity Pfj. Over the frequency interval WA/ it will therefore have 
a total power of PjvVVA/. We may then calculate the source (i.e. 
noise-subtracted) rms amplitude as 

& = (S -P^,W^fY^\ (A7) 

(assuming the noise level P^ is well determined so we can neglect 
any uncertainty in this). Applying the usual rules for transformation 
of variances we find the standard deviation of the estimator is 



sd[a-] 



1 

2& 



M ^ 



1/2 



(A8) 



Of course, we do not know the values of Pj, only their estimates 
yj, so we obtain a reasonable estimate of the uncertainty on the rms 
with the formula 



sd[(7-] 



A/ 



2(tM'/2 



1/2 



(A9) 



This is essentially just the "propagation" (by adding in quadrature) 
of the standard deviations of each periodogram value, which are 
themselves ~ y'., through the averaging over M segments and sum- 
mation over W frequencies (each of which is independently dis- 
tributed). 

This is different from the prescription used bv lOleissner et al.l 

( |2004|) . Their formula is only valid when the power is the same 
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